/******* ********************************, 

************************** Appzuuk 

***********************,, Least Square lattice 

************************ Noise Cancelling ************************,«" 

/* Example for ratiometric approach to noise cancelling */ 7 
/define LAMBDA 0.95 



^*********************** OHUHiH ^ 4 
************* ***************** 
***********************« 4AAAi 



void OxiLSL_NC( int 
int 
int 
. int 
int 
int 
int 



reset, 

passes, 

*signal_l, 

*signal_2, 

*signal_3, 

*target_i, 

*target_2) 



{ 



int 

static int 
static float 



i, ii, k, m, n, contraction; 
*s_a, *s_b, *s_c, *out_a. 
Delta_sqr, scale, noisejref; 



*out_c; 



if( reset 

£T out_a > 

out_c 
N 1 factor 
Li scale 



« TRUE) { 
signal_l ; 
» signal_2, 
= signal_3; 
» target_l ( 
= target_2, 
- 1.5; 

» 1.0 /4160.0; 



noise canceller initialization at time t=0 */ 

£ ncfOJ.berr = o.O; 

s ncfoj. Gamma = i # o; 

Q 

Q for(m=0; m<NC CELLS ; m++) { 



ncfmj.err a 
ncfmj .err~b 
ncfmj .Roh_a 
ncfmj .Roh~c 
ncfmj .Delta 
ncfmj .Fswsqr 
ncfmj . Bsvsqr 



0.0; 
0.0; 
0.0; 
0.0; 
0.0; 

0.00001; 
0.00001; 



/*«■ 

for(k«0; k<passes; k++) { 



END INITIALIZATION 



contraction «= FALSE; 
for(m=0; m< NC_CELLS; m++) { 

ncfmj.berrl « ncfmj.berr; 

ncfmJ.Bswsqrl « ncfm] .Bswsqr; 



/* Update delay elements 



noisejref = factor * logfl.o - (*s_a) * scale) 

- logfl.O - (*s_b) * scale) ; 

nc[0].err_a « log(l.O - (*s_b) * scale); 

ncf0].errj> *= log(i.o - (*s c) * scale); 



# 



Ms 



++s_a ; 
++s_b; 
++s_c; 

nc[0] . ferr 
nc[0] .berr 
nc[0] .Fswsqr 
nc(0] .Bswsqr 



noise_ref ; 
noise_ref ; 

LAMBDA * nc{0]. Fswsqr + noise_ref * noise_ref; 
nc(OJ .Fswsqr; 



/* Order Update */ 

for(n°l;( n < NC_CELLS) it (contraction — FALSE); n++) { 

/* Adaptive Lattice Section */ 

m « n-1; 
ii- n-l; 



nc(m]. Delta *« 
nc(m]. Delta +«= 
Delta_sqr « 



LAMBDA ; 

nc[m].berrl * nc[m].ferr / nc[m] •Gamma ; 
nc(m]. Delta * nc[m). Delta; 



nc[n] . f ref 
nc(n] .bref 

nc[n] .ferr 
nc(n] .berr 

nc[n] .Fswsqr 
nc[n] .Bswsqr 



= -nc(m]. Delta / nc[m] .Bswsqrl; 
= -nc[m). Delta / nc[m) .Fswsqr; 

« nc[m].ferr + nc[n].fref * nc(m].berrl; 
= nc[m].berrl + nc(n].bref * nc[m].ferr; 

nc[m]\Fswsqr - Delta_sqr / nc[m] . Bswsqrl; 
nc[m]. Bswsqrl - Delta_sqr / nc[m] .Fswsqr; 



iff (nc[n] .Fswsqr + nc[n] .Bswsqr) > 0.00001 | | (n < 5) ) { 

nc(n). Gamma = nc[m]. Gamma - nc[m].berrl * nc[m].berrl / nc[m] . Bswsqrl; 
if (nc[n) .Gamma < 0.05) nc[n]. Gamma « 0.05; 
if (nc(n] .Gamma > 1.00) nc[n). Gamma = 1.00; 

/* Joint Process Estimation Section */ 

nc[m].Roh_a *= LAMBDA; 

nc(m).Roh_a += nc(m).berr * nc(m].err_a / nc(m). Gamma ; 
nc[m].k_a « nc(mJ.Roh_a / nc(m) .Bswsqr; 
nc[n).err_a « nc[m].err_a - nc[m].k_a * nc[m].berr; 

nc(m].Roh_c *«= LAMBDA; 

nc(m].Roh_c +« nc[mj.berr * nc[m].err_b / nc[m]. Gamma ; 

nc[m).k c « nc[in).Roh_c / nc[m] .Bswsqr; 

nc[n].err_b « nc(m].err_b - nc[m).k_c * nc[m].berr; 

} 

else { 

contraction «= TRUE; 
for(i=n; i<NC CELLS; i++) { 



nc(i; 


. err_a 




0.0; 


nc[i' 


. Roh_a 




0.0; 


nc[i; 


. err_b 




0.0; 


nc( i; 


.Roh c 




0.0; 


nc[ i; 


.Delta 




0.0; 


nc[ i' 


. Fswsqr 


ss 


0.00001; 


nc(i; 


. Bswsqr 


tB 


o.ooooi; 


nc( i; 


. Bswsqrl 




o.ooooi; 



• # 



} 

) 

} 

*out_a++ = (int) ( (-exp(nc[ii] .err_a) +1.0) / scale) ; 
*out_c++ •= (int) ( (-exp(nc(ii] .err_b) +1.0) / scale) ; 

) 

) 

/********************** Least Square Lattice ****** 



